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BACKGROUND OF THE INVENTION 



I, Field of the Invention 

10 TECHNICAL FIELD 

. The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
.communications with data rates up to multiple Tl (1.544 Mbps) and 
higher {>100 Mbps), and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
means for a new approach to the design of waveforms and filters 
using mathematical formulations which generalize the Wavelet 
concept to communications and radar. 
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BACKGROUND ART 



II, Description of the Related Art 



Multi-resolution waveforms used for signaling and/or 
5 filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency^ with scale 
and shift properties of multi-resolution over the time- frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
10 Wavelet as special cases. The emphasis will be on digital 
design and applications with the understanding that these multi- 
resolution waveforms are equally applicable to analog design and 
applications. 

15 Background art consists of the collection of waveform and 

filtering design techniques which can be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 

20 filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and primarily theorectical studies, C3) combinations of CI and 
C2 for greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques to control the 

25 noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 
scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 

30 techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 
solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include adaptive 

35 equalization for communications, adaptive antenna filters. 
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Wiener filters, Kalman filters, and stochastic optimization 
filters • 

Category CI common examples of LS digital filter design are 
5 the eigenvalue algorith m in ^^A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" [7] and 
the Remez-Exchange algorithm in ^'A Computer Program for Designing 
Optimum FIR Linear Phase Filters^ . - f8^ — The eigenvalue algorithm 
is a direct LS minimization and the Remez-Exchange can be 

10 reformulated as an equivalent LS gradient problem through proper 
choice of the cost function. Both LS algorithms use the FIR 
(finite impulse response) d igital samples as the set of design 
coordinates. LS design metrics are the error residuals in 
meeting their passband and stopband ideal performance as shown 

15 in FIG. 3, Category C2 common examples of analytical waveforms 
and filters for system applications are the analog Chebyshev, 
Elliptic, Butterworth, and the digital raised-cosine, and square- 
root raised-cosine. For theorectical studies, common examples 
are polyphase multirate filters, quadrature mirror filters 

20 (QMF) , and perfect reconstruction filters. Although these 
theorectical studies have yet to yield realizable useful filters 
for system applications, their importance for this invention lies 
in their identification and application of ideal performance 
metrics for filter designs. Category C3 common digital example; 

25 is to start with the derivation of a Remez-Exchange finite 
impulse response (FIR) filter and then up-sample and filter with 
another bandwidth limiting filter. This results in an FIR over 
the desired frequency band that is larger than available with the 
Remez-Exchange algorithm and with sidelobes that now drop off 

30 with frequency compared to the flat sidelobes of the original 
Remez-Exchange FIR. A category C4 common example is to select 
the free parameters of the category C3 filter in order to 
minimize the signal to noise power ratio of the data symbol (SNR) 
losses from the ISI, ACI, and the non-ideal demodulation. A 

35 second common example is start with a truncated pulse whose 



3 



length is short enough compared to the symbol repetition 
interval, to accommodate the timing offsets without significant 
impact on the ISI and ACI SNR losses. This shortened pulse can 
then be shaped in the frequency domain, 

5 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
Gaussian minimiom shift keying (GMSK) waveform. This is a 

10 constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating HPA. 
Category C7 common examples are the adaptive equalization 
filter for communication channels, the adaptive antenna filter, 
and the Kalman filter for applications including target tracking 

15 and prediction as well as for equalization and adaptive antennas. 

Minimizing excess bandwith in the waveform and filter 
design is a key goal in the application of the C1,...,C6 design 
techniques for communications and radar. Excess bandwith is 

20 identified as the symbol a in the bandwidth-time product 
BTs=l+a where the two-sided available frequency band is B and the 
symbol repetition interval is Tg. Current performance capability 
is represented by the use of the square-root raised cosine (sq-rt 
rc) waveform with a= 0.22 to 0.4 as shown in FIG. 46. The goal 

25 is to design a waveform with a=0 within the performance 
constraints of ISI, ACI, passband, sideband, and passband 
ripple. This goal of eliminating the excess bandwidth corresponds 
to the symbol rate equal to the available frequency band l/Tg = 
B/(l+a) == B for a=0. This symbol rate l/Ts=B is well known to be 

30 the maximum possible rate for which orthogonality between 
symbols is maintained. A fundamental performance characteristic 
of our new waveform designs is the ability to eliminate excess 
bandwidth for many applications. 
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Scope of invention will include all of the waveform and 
filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 
will be on the category C5 Wavelets to establish the background 
5 art since Wavelets are multi^resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 
sideband, and passband ripple. Category C2 theorectical studies 

10 also eliminate the excess bandwidth. However, they are not 

multi-resolution waveforms and do not have realizable FIR design 
algorithms. So the emphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 

15 Wavelet background art relevant to this invention consists 

of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
frequency (t-f) signal space, and the design of Wavelets using 
the iterated filter construction. Wavelets are waveforms of 

20 finite extent in time (t) and frequency (f) over the t-f space, 
with multi'-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
defined by equations (1) and (2) as per Daubechies^s ^^Ten 
Lectures on Wavelets^^, Philadelphia :SIAM, 1992 

25 roforoncc — [1] ■ 

Continuous wavelet 
Discrete wavelet 

30 WaM = \i""wi—) (2) 

a 

where the two index parameters "a,b" are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 
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The y is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

indexed Wavelet \|;ab. The scale factor "'|a|"'^^" has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
5 "a,b". Norm is the square root of the energy of the Wavelet 
response. The Wavelets \t/a,b and \\f are localized functions in the 
t^f space which means that both their time and frequency lengths 
are bounded. The discrete Wavelet has the time "t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
10 uniformly sampled at '^T" second intervals. 

Wavelets in digital t-f space have an orthogonal basis that 
is obtained by restricting the choice of the parameters ""a^b" to 
the values a-2"^ , b=qM2P where 'p,q" are the new scale and 
15 translation parameters and "M" is the spacing or repetition 
interval T3=MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale '"p". Wavelets at "p,q" 
are related to the mother Wavelet by the equation -[4+ 

20 "^'^ J r\ ^ J ^3 J 

-2-^^«^(2-Pn-qM) 



where the mother Wavelet is a real and even function of the 
sample coordinates. The orthonormality property means that these 
25 Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1". 

E^'p.qM'k.m =1 iflfboth p = kandq = m 

n 

= 0 otherwise (4) 

Wavelet representation of a digital t-f space starts with 
30 selecting an N sample time window of a uniform stream of 
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digital samples at the rate of 1/T Hz (1/second) equivalent to a 
"T" second sampling interval. The N point or sample t-f space 
in FIG, 1 illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
5 construction. 

The t"f space in FIG. 1 is partitioned or covered or tiled 

by a set of Wavelet subspaces ^^,/? = 0,l,,..,m-l} 1 where N=2^ Each 

Wavelet subspace Wp 1 at scale '^p" consists of the set of Wavelet 
10 time translations {q = 0,1,..,, N/2^^^-1} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 
each subspace are mutually orthogonal with respect to the time 
translates. This N-point t-f space extends over the time 
interval from 0 to (N-DT 2 where T is the digital sampling 
15 interval, and over the frequency interval from 0 to (N-1) 3 in 
units of the normalized frequency fNT 4. 

The iterated filter bank in FIG, 2 is used to generate the 
Wavelets which cover the t-f space in FIG. 1. Each filter stage 
20 5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . . Output 6 of the LPF is subsampled by 2 which is 
equivalent to decimation by 2. This t-f space space is an N- 
dimensional complex vector metric space F 7. At stage m in the 
iterated filter bank, the remaining t-f space F^.j 8 is 

25 partitioned into V^^j 9 and the Wavelet subspace Wm+i 10 . 

Scaling functions and Wavelets at each stage of this filter 
bank satisfy the following equations 

30 

<p(«)=2-^^X^^(2«-g) Vq 

1 
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where <p is the scaling function, V}/ is the Wavelet, HPFp 

coefficients are {h^^Vq}, LPFp coefficients are { g^, Vq } , and 

the equations apply to the stages 0, 1, ...,m-l . . Identifying the 
5 scale parameter and using the previous Wavelet formulations 
enable these equations to be rewritten for stages p=0, 1, ...,m-l 
as 

q 

10 

For our application the HPFp and LPFp are quadrature mirror 
filters (QMF) with perfect reconstruction. This means they cover 
the subspace Vp with flat responses over the subband frequency 
including the edges of the frequency subband, and the HPFp 
15 coefficients are the frequency translated coefficients for the 
LPFp: {gq-(-l)% , Vq}. 

Wavelet design using iterated filter bank starts with the 
selection of the scaling functions. Starting with a primative 
20 scaling function such as the one proposed by Daubechies — [1] , 
J_one can use the iterated filter construction given by equations 
(5) and (6) to derive successive approximations to a desired 
scaling function which has properties that have been designed 
into it by the selection of the filter coefficients {g^,Vq} at 

25 each level of iteration. The Wavelets can be derived from these 
scaling functions using the iterated filter construction or 
scaling equations (5) and (6) . 

Another use of the iterated filter construction is to 
30 design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution analysis and 
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synthesis as illustrated by Coifman's Wavelets in ^^Wavelet 
analysis and signal processings^ > — Dcoign CKamploo abound in 
^fefee — opon — literature — and the — rocont publicationo — — [4-J-7 — [5] ^ 
— arc roprcocntativo illuotrationa of the divoroo applicationo 
5 of Wavclcto. 

SUMMARY OF INVENTION 

10 SUMMARY OF THE INVENTION 

This invention is a new and novel design of Wavelet multi- 
resolution waveforms and filters in the frequency domain with a 
property which provides a single waveform design for all of the 

15 waveforms at multiple resolutions and scales. Wavelets are 
designed to meet the specific application requirements, are 
complex, are defined by frequency harmonics, include a frequency 
translation parameter in addition to the scale and translation 
parameters in current Wavelets, and have the property that the 

20 set of design coefficients in the frequency domain remain 
invariant for all possible multi-resolution scales, and time- 
frequency translations. Frequency design harmonics and the 
frequency translation capability enable the waveforms at multiple 
resolution scales to be derived from the single waveform design 

25 by scaling the dilation, translation, and frequency translation 
parameters using the set of frequency design harmonics. Wavelet 
waveform design is illustrated by Matlab 5.0 code to design a 
linear filter waveform using an iterative least-squares 
eigenvalue approach to minimize the non-linear least-squares cost 

30 function, and to scale this waveform design for multi-resolution 
application specified by the dilation, time translation, and 
frequency translation parameters. Additional results are given 
for a constant amplitude minimum-shift-keying bandwidth-efficient 
modulation waveform and for a synthetic aperture radar waveform. 
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This invention is a new approach to the design of multi- 
resolution waveforms that improves their performance for 
engineering and scientific applications, and in particular for 
5 applications to communications and radar. Current practice is to 
1) design traditional multi-resolution waveforms in the time 
domain using metrics which specify frequency performance and 
without consideration of Wavelet propertiesT — aa d or to 2) design 
Wavelet multi-resolution waveforms using the time-domain iterated 
10 multi-resolution filtering approach with a set of scaling 
functions used to perform the filtering and without direct 
considerations of the frequency performance. These two separate 
design approaches yield fundamentally different waveforms. Our 
new invention provides a means to combine these two approaches 
15 to generate a new multi-resolution waveform with the best 
properties of the traditional multi-resolution waveforms and the 
Wavelet multi-resolution waveform. 

Our invention introducca tho new innovationo for tho dooign 
— ettae — aew — waveform — that — 1-) — include — the — tS-^T — — aftd — QMF 
rcquircmonto — into — the — dooign — algorithm — along — with — non linear 
modifications to tho traditional paooband and otopband froquoncy 
roquiromonto which arc used in the Rcmcz - Exchango and oigonvaluc 

LS algorithmo in rcfcrcncGO — W-j — — roopcctivoly/ ^ — ocloct 

^fefee — frequency harmonics — ete — the — dooign coordinatoa / — -3-) — u-se — aft 
appropriate — ouboot — e# — the — available — Fourior — domain — froquoncy 
harmonico ao the dooign coordinatoo v^yith the property that thio 
ouboot — e^ — coordinatoo — i-s — a — baoio — fe^ — the — multi - reoolution 
waveform dooign, and 4) — dooign the multi - rooolution vmvoform for 

no cxccoo bandwidth a-Q. 

Thooo — innovationo — support — the — development — of our — aew 

waveformo — h) that — a^?e — goncraliaationo . of — Wavoloto — — the 

froquoncy — domain — uoing — a — mcano — which — makoo — them — uocful — fe^? 
communicationo — and radar — ovor their t - f — opacoy — 3-) — that — can bo 

10 
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dooignod — fee — fee — aft — orthonormal — baoi£3 — e* — orthonormal — sefe — e# 
coordinatGO — over — febe — available — f roqucncy — bandwidth — with — fefee 

implicit — property — that — tfee — oxoooa — bandwidth — vaniohoo a-0>> — fee 

vfithin the acGuracico allowed by communication and radar dcoign 
5 implomcntationo , — and 3) — that have multi - oealc propertieo whieh 
allow a oinglc do waveform dcoign to be uocd to uniquely dofinc 
the complete — oct of new vraveformo over the t-f opaco for multi - 
ocalc applications Specific dcoign algorithm cxampleo developed 
in thio — invention dioclooure arc the LS dcoign algorithmo — that 
10 uac eigenvalue and gradient — oearch tcchniquco — in the — frequency 
domain — fee — find — fefee — boot — wavcfom — dcoign — which — minimi zeo — fetie 
corresponding LS error rcoidual coot function that io a weighted 
linear oum of the rcoidual error mctrico. 

15 Waveform performance calculated ift thio invcnoion 

dioclooure from the application of these LS deoign algorithmo, 

dcmonotratco the capability of thio invention to provide a means 
fee — deoign — wavcformo — that — et^ — improved — over — current — practice . 
There — aa?e — other — algorithmic — design — algorithms — which — esrR — fee 

20 realized by invention. Aft — example — given — ift — thio — invention 

disclosure — is — fefee — modification — of — theoo — LS — algorithmo — fe*^ 

application — to conotant — amplitude — bandwidth — efficient (BEM) 

wavcformo for communications. 

25 BRIEF DESCRIPTION OF THE DRAWINGS AND 

THE PERFORMANCE DATA 

BRIEF DE S CRIPTION OF DRAWINGS MID PBRFOIiMANCE DATA 

The above-mentioned and other features, objects, design 
30 algorithms, and performance advantages of the present invention 
will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
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data wherein like reference characters and numerals denote like 
elements, and in which: 

FIG* 1 is a Wavelet eRr-N-point t-f space extending over the 
5 time interval (0, (N-l)T] and the frequency interval (0, (N- 
1)/NT], which is tiled or covered by a set of orthonormal 
Wavelets at the scales p==0, l,...,m-l . 

FIG* 2 is a Wavelet iterated filter bank used to generate 
10 the set of Wavelets which tile or cover the t-f space in FIG. 1. 

FIG. — 3 illuQtratoo the power opcctral density — (PSD) — of our 

now waveform at do, and the oct of otopband and paooband design 

rcquircmcnto . 

15 

FIG. 3 is the power spectral density (PSD) template for the 
the stopband and passband used to construct the LS error metrics, 

20 FIG. 4 is the flow diagram of the LS metrics and cost 

functions and the final cost function used to find the optimal LS 
solution for our the Wavelet waveform. 

FIG. — 5 — ia — flow — diagram — fea? — the — hS — rccuroivc — oolution 
25 algorithm to find the — optimal harmonic coordinatoo — that dofino 
the LS oolution for our waveform. 

FIG. 5 is the Matlab 5.0 code for the LS algorithm used to 
design the Wavelet waveform in FIG. 6. 

30 

¥^LG-. — 6 ploto the do PSD in dB unito for tho wavoformo vo . 
the normalized frequoncy fT^ for the — vjavcform dooigncd with the 
— algorithms — aad — #ea? — the — oq - rt — e — waveforms — with — excooo 
bandwidth parameter — 0.22, — 0.40. 
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FIG. 6 plots the dc PSD in dB units for the frequency 
response of the mother Wavelet at dc and the square-root raised- 
cosine (sq-rt r-c) waveforms with excess bandwidth parameter a= 
5 0.22, 0.40 vs. the normalized frequency in units of symbol rate 
1/Ts. 

FIG. — 7 plots tho dc PSD in dB unito for our wavoformo vo. 
the normaligcd frequency fT^ for tho now wavoform dooignod with 
10 modif icationo — to tho LS algorithms — for application to — constant 

amplitude BEM communicationo, and for a Gauooian minimum shift 

keying — (CMSK) wavoform. 

FIG. 7 plots the dc PSD in dB units for the frequency 
15 response of the mother Wavelet at dc and an optimized Gaussian 
minimum shift keying (GMSK) waveform vs. the normalized 
frequency fT s. 

FIG. 8 plots the amplitude of the dc radar ambiguity 
20 function vs. the normalized frequency fTp and normalized time tTc 
where the pulse time Tp and chip time Tc are both equal to the 
communications symbol interval Tg for this example, for the new 
waveform and an unweighted chirp waveform. 

25 
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DI S CLOSURE OF IMVEaJTION 

DISCLOSURE OF THE INVENTION 

35 



13 



Wavelet multi-resolution waveforms incorporate a frequency 
translation, are complex, are designed in frequency harmonic 
coordinates, a single mother Wavelet design is used for all 
scales, and in the t-f space are realized as FIR filters and 
5 waveforms . 

lie wThese Wavelet waveforms in this invention disclosure are 
generalizations of Wavelets in t-f space which enable them to be 
useful for communications and radar applications. This 

10 generalization is accomplished by 1) the addition introduction 
of a frequency translation, by 2) relaxing — changing the 
orthonormality condition in equation (4) to apply to waveforms 
within the same space {q} and over the scales {p} with the 
inclusion of the frequency translation, and 3) by their 

15 characterization and design in the Fourier domain. With frequency 
translation the analytical formulation of these new waveforms as 
a function of the baseband or mother waveform centered at dc( dc 
refers to the origin f=0 of the frequency space) becomes 

20 

New waveform ao a function of dc waveform 

yfp^,^{n) = 2-P^K2'^n-qM)e*'^^^'^^ (7) 

25 where fc(P/r) is the center frequency of the frequency translated 
dc waveform, at scale "p" and frequency index "r". The purpose 
of the frequency index ^'r" is to identify the center frequencies 
of the waveforms at the scale ^'p" in the t-f space. The dc or 
baseband or mother waveform is the generalization of the mother 

30 Wavelet for multi-resolution waveforms - 

These ' waveforms satisfy the complex orthonormality 
equations 
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Orthogonality oquationo 

V*iL =1 iflf p = kandq = m and r = V 

= 0 otherwise (8) 



where is conjugation, and are generalizations of the 

5 orthonormality equations for the analytical Wavelets in (4) . 

— ftew — wQvcf ormo — — equation — W expand — =fefee — Wavolct 

analytical formulation te include — a frequency — variable , 

Wavcleto — are — functiono of the — scale and tranolation paramctero 

10 ^^P/q^^. The concept of an additional parameter to provide an added 
degree of flexibility in their tabulation was f irot introduced by 

^tbe — Coif man — library — e# Wavcleto — — which — uee — aa — oocillaton 

parameter — §^ — tabulation — that — roughly — corrcopondo — fee — fehe 
frequency of oscillation of the Wavcleto. The frequency variable 

15 together with the Fourier domain design — are entirely new meano 
-§e^^ — deriving — these — new — wave f ormo — as — generalization — e^ — feiie 
traditional — Wavcleto — aftd — their — modification — with — the — added 
parameter in [2] . 

20 New These waveforms are generalizations of Wavelets in the 

frequency domain using a means which makes them useful for 
communications and radar in the t-f space. The basis vectors for 
this metric space V consist of a subset of the admissible set 
of scaled and translated waveforms {^p^q^r^^V-^^^^ derived from the 

25 dc waveform Xj/ as per equation (7) . An admissible waveform is 

any . combination that covers F=t-f space. We are interested in 
the Fourier domain representation of the dc waveform ^ in F , 
and in particular in a subset of the discrete Fourier transform 
(DFT) harmonic coefficients over the Fourier domain which we 
30 intend to use as the design coordinates- Starting with the z- 
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transform and continuous Fourier transform, the DFT harmonic 
coefficients are defined by the following equations 

DFT Harmonic Coefficients (9) 
y/{z) = ^\f/{xi) 7^ z-transform 

n 

y/(a)) = ^^(n) e"'"*" Fourier transform 

n 

y^k = ^4^^^^^ DFT harmonic coefficients for \fk 

n 

where 

iff{a>) = (1/ N*) Xe'^'^"^ 

k n 

- sin((^/2- ;;*/N')iV')/ Nsin(ty/2 -7^) 
k 

- 2 [Harmonic interpolation for "k"] 

k 

{k} = DFT frequency or harmonic coefBcients such that 4 - k/N'T 
where f,, is the harmonic frequeuQ^ corresponding to "k** 
= length of y/{ti) 

where ^(fi?) = is the power spectral density of the dc 

waveform. These equations define the frequency representation of 
our new waveforms in terms of the available set of harmonic 
coefficients, which set is considerable larger than required for 
most applications. 

New — Wavelet wave forms ei?e — an — orthonormal — baoio — with — ne 

cxGQOo bandwidth, aro propcrtioo that v<rc will dcmonotrato have . 

Orthomality orthomality and no excess bandwidth arc p roperties 
which are asymptoticly approached by our new wavefoms to within 
design accuracies inherent in communications and radar • To 
proGQod — demonstrate these properties w e need to identify the 
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structure of the dc waveform in V . We start with definitions 
for the parameters and coordinates in the following equation 
(10) . The waveforms derived from the dc waveform will be 
designed to be orthogonal over both time translates ^^MT" and 
5 frequency translates ^'l/ij?MT" which respectively correspond to 
the Wavelet symbol spacing Ts=MT and the adjacent channel 
spacing l/Tg. This means the orthogonal spacing of the waveforms 
in V are at the time^f requency increments (MT, = (Ts- 

1/Ts) . In the interests of constructing our ftew — orthonormal 
10 multi-resolution waveforms to cover V it will be convenient to 
assume that M,L are powers of 2. We need the following 
definitions for the parameters and coordinates. 
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Parameters and Coordinates 



N' = Length of \(i which is an even function about the center and 
which spans an odd number of points or sannples 
= ML + 1 where M, L are assumed to be even fimctions for 

convenience of this analysis 
= Number of points of \|/ 
M = Sanq)ling interval for y/ 

= Spacing of ^ for orthogonality 
L = Length of v|f in units of the sample interval M 

= Stretching of \|/ over L sanq)le intervals 
n = Ho + HjM 

= partitioning into an index n^ over the sanq>le 
length no =0,1,..,,M -land an index n, =0,1...,L-1 
over the sanfq)le intervals 
k =ko -f k,L 

= partitioning into an index k^j over the harmonic 
frequencies k^ = 0,1,...,L-1 corresponding to the 
stretching and an index k^ = 0,1,...,M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for I// 

The harmonic design coordinates are selected using the 
5 following observation. For most applications and in the following 
development it is assximed that the waveforms are spectrally 
contained in the frequency interval 1/Ii^ ^T corresponding to the 
frequency spacing • This suggests the harmonic design coordinates 
be restricted to the subset of L harmonics {ko==0, 1, • . . ,L-1} 
10 covering this spacing. These L harmonics correspond to the 
stretching of the mother waveform over the L repetition 
intervals. 

Obviously, for some applications as will be demonstrated 
15 later, the spectral containment is spread out over several 1/LT 
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frequency increments whereupon one must increase the subset of 
design harmonics to possibly 2L, 3L or larger. 



Pursuing this application using L design harmonics, ^ he — the 
5 DFT equations for the dc waveform in (9) when rewritten in terms 

of the L harmonic design coordinates {^^^j^^, Vkp} become: 
DFT equations for dc waveform 

10 

" y^,y^(Q) WjJ°° harmonic design coordinates 

Q 

\{f{n) =(1/N') E'J'k, Wj!j?" new waveform defined in terms (11) 
of the L harmonic design coordinates {^j^^^Vko} 

15 

It will now bo ohown that It will be proven that the use of 
these L harmonic design coordinates is sufficient — our now 
vjavcform — under time translates to be a basis for the 
corresponding subspace of F— Thio — which m eans these waveforms 
20 provide a complete set of coordinates to describe this subspace. 
tfe — The proof w ill use the 

will uoG — the theorems of Karhunen-Loeve and Mercer and will 
limit the demonstration to the dc waveform for simplicity and 
without loss of generality. Wo start by conoidoring C onsider the 

25 expansion of a random complex sequence {z(n), Vn} in a series of 
waveform coordinates consisting of time translates of the mother 
waveform tfr . The sequence {z(n), Vn} is a zero-mean stationary 
random process which is orthonormal over the sample interval '^M" 
and has a frequency spectrxam which is flat and extends over the 

30 frequency range 1/LT which is centered at baseband corresponding 
to a zero frequency. This means the {z(n), Vn} cover the 
subspace of V corresponding to the scale of our new dc waveform 
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^ and its time translates { ^(w-^M) - ^^(w), Vq } . In addition, 

the V is now considered to be extended over a time interval 
which is relatively large compared to the N-dimensional t-f space 
in FIG, —1,2 to avoid end-effects on the analysis. We start by 
5 approximating the sequence {2(n), Vn} by the {2(w),Vn} where: 

2(D)- SZ, Kn-qM) 

4 
4 

where the complex coefficien ts {z^ }are 

derived from the original sequence using the ( 12 ) 

orthogonality properties of i// 
Z, =i:z(n) v^*(n-qM) 

n 

= i:z(n)i/,(ii) 

D . 

10 The following equations prove that the coefficients {Zq, Vq} are 
orthonormal : 



Z^Z,.* = 2:z(An-qM)z*(An-q'MV,(AIl)^^^^. ♦(An) (13) 

Ad 

= S|v^q(n)f since the sequence {z} is orthonormal 

n 

for qM time translates 
= with normalization of the energy of ^ 

15 

Equations (12) and (13) together prove the Karhunen- 
Loeve's theorem which proves the following equation for the 
accuracy in approximating the stochastic sequence {z(n), Vn} by 
20 {z(n),Vn}, This accuracy is expressed by the expected '"E(o)" 
squared error "(o)" in this approximation: 
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(14) 



5 We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 



The right hand side of equation (14) when set equal to zero 
expresses Mercer's theorem so our goal is to prove Mercer's 
theorem. To do this we use the DFT of \f/ in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
15 (14) . We find 



proves that the new waveform coordinates { Vq } are a basis for 



the original sequence {z(n), Vn} which is our goal. 




where T = (l/Lf^Wf'^ 



= sin(;r Ak9)/Isiii(;r Akfl /L) 
= 1 for Ako = 0 
= 0 otherwise 



(15) 



This proves that 



1 - +niM)f = 0 Vno 



20 which proves that equation (3^14) reduces to 



E{z(n)-zinf\= 0 



(16) 



which as per the above proves that the set of multi-resolution 
waveforms is a basis. This proof easily generalizes to the 
multi-resolution waveforms . at all of the scales {p} and time 
5 translates {q}, and to the expansion of the harmonic design 
coordinates over 2L, 3L, as required by the application. 

New — Next we demonstrate that these w aveforms have multi- 
scale properties^^ that wo will dcmonotrato. Wc otart vfith the 

10 oboorvation — that — tbe — The Fourier domain design for the new 
waveforms provides a natural and easy way to derive the complete 
set of waveforms { t^^^,^ , Vp,q } for the space V from the design of 
the dc waveform ^ by .using the same invariant set of Fourier 
domain harmonic design coordinates {%^,Vko} derived for the dc 

15 multi-resolution waveform. This demonstration requires that we 
show 1) how the multi-scale transformations are implemented with 
the design in the Fourier domain, and 2) how the waveform design 
remains invariant under scale changes. 

20 First consider the multi-scale transformation which derives 

the waveforms at the scale and shift parameters "p,q" for the dc 
waveform at scale ^^p=0" and centered at the origin ^"q-0". We 
begin by extending the parameters and coordinates in equation 
(10) , to include both scaling and subsampling or decimation in a 

25 form that is equivalent to the iterated filter bank construction 
which is used to derive the waveform using the filter scaling 
functions. Starting with the coordinates at scale "p=0" the 
parameters and coordinates at scale "p=p" are given by the 
equations: 

30 
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^^£=0 parameters and coordinates 

n = no + ni M (17) 

no = ao + ai 2 + ... + am-i 2^"*^ 

= M points 
M = 2^ 

ni = bo + bi 2 + ... +bi-i 2^"^ 

= L points spaced at M sample intervals 
L = 2 

p-p parameters and coordinates 

2-P n(i2^) - scaled by "2"^" 

and subsampled or decimated by 2^:1 

= no(p) + ni(p)M 

15 2-P no(>^2^) = scaled by ^^2"^" 

and subsampled or decimated by 2^:1 
= no(p) 

= ap + ap.i2 + ... + ap+,a-i 2^^"^"' 

= M points spaced at 2^ sample intervals 

20 2'P ni(>l2^) = scaled by ^^2"^" 

and subsampled or decimated by 2^:1 
= ni(p) 

- bp + bp^i 2 + ... + bp.i-1 2^^^-^ 

= L points spaced at M2P sample intervals 

25 

together with the observation that the sampling interval "T" is 
increased to ^'2^T'' under the scale change from "p-0" to ''p=p" and 
subsampling or decimation from "1:1" to "2^:1", Combining these 
30 equations with the analytical formulation in (7) and the Fourier 
domain representation in (11) enables the waveforms at the 
parameters ''p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 
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1. (n<r.^-n^J^^ '2;r/,(p.rMp)2Pr (18) 



5 

for all admissible scale, translation, and frequency index 
parameters ''p,q,r". 

Next we need to demonstrate that the frequency domain 
10 design in (11) remains invariant for all parameter changes and in 
particular for all scale changes. This multi-scale property 
expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
15 frequency translation indices. This multi-scale invariancy means 
that the design for a M=16 channel filter bank remains the same 
for M=100 or M=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the m^^ — multi- 
20 resolution waveform at scale ^'p" with the other parameters set 
equal to zero for convenience ''q=0, r=0" and without loss of 
generality. The Fourier domain frequency response ^(/) can be 
evaluated starting with the original formulation in equation (7) : 



25 



30 
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DFT at 



"p, q=0, r=0 



(19) 



ko n(p) 

sin(;r (£Nf'2^T-ko)) 
N'sin(;r (fN*2PT.ko)/N')^ 

[Harmonic interpolation for " ko "] 

5 

This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 

harmonic coefficients {^j^^^Vko} and the stretching of the time 

interval to ^^2^1" corresponding to the scale ^'p=p"- The harmonic 
interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 
"f-1/2^" which means the frequency* time" product remains 
invariant with scale changes as per the fundamental property of 
the waveforms. This means the harmonic interpolation functions 
remain invariant with scale change and therefore the frequency 
response remains an invariant. This demonstrates the waveform 
design is an invariant across the waveform scales which means we 
only need a single design for all scales or ' resolutions of 
interest - 

LS design algorithms for new waveform will be described to 
illustrate the advantages our ftew — waveform has over current 
designs. The two LS algorithms described are the eigenvalue and 
25 the gradient search which respectively can be reduced to 

algorithms which are equivalent to ^fefee original — current 

eigenvalue [7] and Remez-exchange [8] w aveform design algorithms 
for application to a uniform filter bank. We consider the t-f 

25 
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space which is spanned by a uniform polyphase filter bank 
consisting of M channels at the frequency spacing fw==l/MT where T 
is the digital sampling interval, and the filter waveform FIR 
time response is stretched over L sampling time intervals Tg. 
5 This polyphase filter bank is ideally decimated which means the 
filter ouput sample rate l/Tg is equal to the channel-to-channel 
spacing Ts=MT, equivalent to stating that there is no excess 
bandwidth a=0. Our design for this topology is immediately 
applicable to an arbitrary set of multi-resolution filters 
10 through the scaling equation (18) which gives the design of our 
waveform at arbitrary scales in terms of our design of the dc 
waveform. 

For this polyphase filter bank used to construct the dc 

15 waveform or filter impulse response, our LS example design 
algorithms will use 5 metrics consisting of the 2 prior art 
passband and stopand metrics, and the 3 new metrics consisting 
of the ISI, ACI, and QMF, and solve the LS minimization problem 
using as design coordinates the subset of harmonic coordinates 

20 which are a basis. Since our S — two example LS design algorithms 
only differ in the use of an eigenvalue LS optimization and the 
use of a gradient search LS optimization, the flow diagrams for 
the construction of the cost functions and the solution for the 
optimal waveform will be identical. However, there are 

25 differences in the construction of the cost functions from the 
respective metrics and in the iterative solution mathematics. 
Both LS solutions for the harmonic design coordinates minimize 
the weighted sum of the error residuals or cost functions from 
the 5 metrics. These design coordinates for are the Fourier 

30 harmonics {i/jt^,Vko} for {ko = 0, 1,..., L-1 } . Resulting algorithms 
are easily extended to the applications requiring the design 
coordinates to cover 2L, 3L, ... harmonics. 
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Frequency domain design coordinates are related to the 
waveform time domain digital samples or coordinates as follows. 



Mappings of timo — ^ — frequency to time (20) 

5 

Time domain design coordinates Vn}are real and 

symmetric and can be represented by the reduced set {ht(n), 
n=0,l,...,ML/2} 

ht(n) = t/r(0) for n=0 
10 = 2ffr(n) for n=l,2,...,ML/2 

= time domain design coordinates 



Frequency domain harmonic design coordinates {^j^^^VTcq} are 

real and symmetric and can be represented by the reduced 
15 set {hf{k), k=0,l, ,...,L-1} 

hf(k) = i/Tf^^ for k=ko=0 

= 2if/f^^ for k=ko=l,2,...,L-l 

= frequency domain design coordinates 

20 Mapping of the frequency coordinates {hf{k), k=0, 1,...,L-1} 

into the time coordinates {ht{n), n=0, l,...,ML/2} is defined 
by the matrix transformation 
ht = B -hf 

where 

25 hf = (hf(0), ...,hf {ML/2) )^ transpose of column vector 

ht = (ht(0), ...,ht (ML/2) ) ^ transpose of column vector 
B = (ML/2 +1) X L matrix 

= [ Bkn ] matrix of row k and col\imn n elements 5^ 
Bkn = 1 /ML for n=l 

30 =2 cos(27ckn/ ML) otherwise 
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Mappingo of the frequency coordinatQa — R^ ^(k) i — k^O, ly .../L - l] 
into — ^he — tifae — coordinatoo — f^^^-f^ii-T — n^O, 1/ .,.,ML/2) — and vioe 
voroa are defined by the matrix tranoformationo 
B — fef 



where 

h ^ = (h^(0), .„^h^(ML/2) )^ — tranoposo of coliimn vector 

(h^(0), „.,h^(ML/2) )^ — tranopooe of column vector 

B — - (ML/2 11) X L matrix 

10 ^ [ B^ tft— j — matrix of row k and column n clemonto Bj tft 

- 1 /ML + 1 '■ for n-l 

= 2 coQ(27tkn/ ML + 1) otherwioe 



• B" ^ - L K (ML/2 11) matrix 

— \rB^i — matrix of row k and column n clcmento -B^ 

15 9^-^ for n-1 

- 2 coo (27Ckn/ MLH) otherwioe 



and 



-B" ^ - I L X L identify matrix 

-B"B" ^ - I (ML/2 + 1) X (ML/2 i l) identity matrix 



20 

wherein the N^=ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
ease of implementation for the example Wavelet being considered 
with a sample at the mid-point that makes the number of samples 
25 N' an odd number. 



Passband and stopband metrics and cost functions are 
derived with the aid of FIG. 3 which defines the power spectral 
30 density (PSD) parameters of interest for the passband and 
stopband of the PSD ^((d) for communications applications. 
Requirements for radar applications include these listed for 
communications. Referring to FIG. 3 the passband 11 of the 
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waveform PSD is centered at dc {f=0) since we are designing the 
dc or baseband waveforin, and extends over the frequency range % 
extending from -cOp/2 to +o)p/2 12 in units of the radian 
frequency variable (0=27cfT 13 where T is the digital sampling 
5 interval defined in FIG. 1. The frequency space extends over the 
range of f=-l/2T to f=+l/2T which is the frequency range in 
FIG 1 translated by -1/2T so that the dc waveform is at the 
center of the frequency band. Quality of the PSD over the 
passband is expressed by the passband ripple 14. Stopband 15 

10 starts at the edge 16 of the passbands of the adjacent channels 
+ /-0)a/2 16 and extends to the edge of the frequency band (0=+/- 
n 17 respectively. Stopband attenuation 18 at +/-(0a/2 
measures the PSD isolation between the edge of the passband for 
the dc waveform and the start of the passband for the adjacent 

15 channels centered at +/-a)s 19- Rolloff 20 of the stopband is 
required to mitigate the spillover of the channels other than the 
adjacent channels, onto the dc channel, Deadband or transition 
band 21 is the interval between the passbands of contiguous 
channels, and is illustrated in FIG. 3 by the interval from a)p/2 

20 to a)a/2 between the dc channel and adjacent channel at (Oa* 
Waveform sample rate cOg 22 is the waveform repetition rate. 
For the LS example algorithms, the waveform sample rate is equal 
to the channel-to-channel spacing for zero excess bandwidth. 
Therefore, 1/Ts= <x)s/2nT-l/MT which, can be solved to give C£)s=27c/M 

25 for the radian frequency sampling rate of the filter bank which 
is identical to the waveform repetition rate. 

We start by rewriting the DFT equations for the dc waveform 
in (11) as a function of the {hf{k), k=0, 1, L-1} 

30 
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y/{o)) = ^^(n)cos(nfi>) 

n 

= c* B hf using (21) and the definition of the vector "c" (21) 
c = ( 1 * cos(fi>), cos( (ML/2) Q)) y transpose of column vector 

which is equivalent to the equation for yfico) in (9) expressed as 
5 a linear function of the {hf(k), k=0, 1,...,L-1} . However, — this 
. functional — form ouggcotod by the eigenvalue — formulation — — ats- 

more — convenient — fee — analyse . An ideal ^^c" vector ^^c^" will be 

introduced for the passband and the stopband in FIG. 3y in order 
to identify the error residual S\f/(Q)) at the frequency ^^co" in 
10 meeting the ideal passband and stopband requirements. The ideal 
PSD is flat and equal to "1" for the passband, and equal to ^'0" 
for the stopband. We find 

15 Error residuals for passband and stopband (22) 

Cr = (1,1,...,!)^ passband ideal ^^c" 
= (0,0,...,0)^ stopband ideal ^^c" 
Sc = Cr - c error vector 

- Residxial error in meeting the ideal spectrum at "c?" 

The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
25 (or equivalently, the LS algorithm) design algorithms— 
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Passband and stopband metrics 



J (band) = — — j|<yt/(fi;)|'' da) Eigenvalue 



band 

= Ikv' II 



Eigenvalue 
LS 



(23) — continued 



-where 



10 



15 



band 



band = [0, (Op ) passband 
= (©s , 71] stopband 

— ((B*&&*B)da> 

band 

= L X L matrix 
= ( 8\|/(G)i) 5\|/(o>u) )^ 
= vector of error residuals at the 

frequencies ©i , (Ou across the band 

norm or length of the vector (o) and which 

includes a cost function for the errors of 
the individual components 



R 



M - 



20 



25 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
coordinates {hf{k), k==0, 1,...,L-1} whereas with the LS approach it 
is sufficient to give the LS metrics as vector norms with 
imbedded cost functions or an equivalent formulation. 
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QMF metrics express the requirements on the deadband that 
the PSD' s from the contiguous channels in FIG. 3 add to unity 
across the deadband [a>p , cOs] in order that the filters be QMF 
filters. By suitable modification of the error vector 5c, the 
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previous construction of the passband and stopband metrics can be 
modified to apply to the deadband. Thio io a Wavolot roquiromont 
on the — f rcquoncy coverage — of the Wavelet baois — as — obocrvcd in 
Tabic 1. ^We find 



5 
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Deadband metrics (24) 



J (deadband) = (I^V^(^)f Eigenvalue 



deadband ^^^^^.^^ 

(24) continued 
= /i^R hf Eigenvalue 

= )Syf II LS 

15 where 

8c = Cr - c(a)) - C(7C/M - (0) 



where c(<o)=c as defined in (22) and (23), and c(7t/M - co) = c at 
20 the offset frequency ""n/M - ©" corresponding to the overlap of 
the contiguous filters over the deadband. 

Orthonormality metrics measure how close we are able to 
designing the set of waveforms to be orthonormal over the t-f 

25 space, with the closeness given by the ISI and the ACI. ISI io 
^fciie — non - orthogonality — error — between — channel — output — oamploo 

ocparatcd by multiplco — of the oampling interval 1/MT occondo 

where T io the oamplc time and M io the interval of contiguouo 
aamplco. — ftGi — i-s — t^ie — non -orthogonality — error — bctvfccn between 

30 channel — output — oamplco — within — a — channel — aed — fcfee — oampleo — ift 
adjacent — channolo — at the — oamo oample time and at oamplc timeo 
separated by multiplco — of the oamplc interval. Ao obocrvcd ao 
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noioc contributiono within each oamplc in a given channel, the 

ISI io the noioc contribution due to the other received wavcformo 
at the different timing offoeto correoponding to multipleo of the 

oampling interval. Likewioc/ the ACI io the noioc contribution 

5 dtte — te — tfee — other wavcformo — iR — adjacent — channelo — at — the — oame 
oampling time and at multipleo of the oampling interval. 
ISI and ACI errors are fundamentally caused by different 
mechanisms and therefore have separate metrics and weights to 
specify their relative importance to the overall sum of the LS 

10 metrics. ISI is a measure of the non-orthogonality between the 
stream of waveforms within a channel as per the construction in 
FIG. 3. On the other hand, ACI is a measure of the non- 
orthogonality between the waveform within a channel and the other 
waveforms in adjacent channels. This means the stopband 

15 performance metric has a significant impact on the ACI due to the 
sharp rolloff in frequency of the adjacent channel, and the ACI 
metric is then a measure of the residual non-orthogonality due to 
the inability of the stopband rolloff in frequency from 
completely eliminating the ACI errors. 

20 

We assiime that the received waveform is identical to the 
filter waveforms and is transmitted at the filter output sample 
intervals equal to i/MT seconds. The second assumption means we 
are assuming the receiver is synchronized with the received 
25 signal. Since there is no information lost by sampling 
asynchronously with the received waveform, we are free to make 
this synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 

30 ISI metrics (25) 

Mapping of hf into \\f 

Vj/ = {v|/{-ML/2), ... , H/(ML/2))^ transpose of column vector 
^ H hf 
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H = (ML+l) X L matrix of elements 

Hjcn 1 for n=l 

= 0.5 cos (27tkn/ML+l) otherwise 
Offset matrix A 

5 A = L x (ML+l) matrix of elements Akn 

Akn =0 for k=l 

= [0 0 \|/(-ML/2) V|/{ {L-k)M+l) ] (25) continued 

12 kM LM+1 

10 

ISI error vector 8E 
5E = L x 1 column vector 
= A H hf 

15 ISI metric 

J(ISI) = 5e^ 8E Eigenvalue 

= Non-linear quadratic function of hf 
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ACI metrics are derived using the ISI metric equations with the 
following modifications. 

ACI metrics (26) 



Mapping of hf into i|/ is the same as developed for ISI 
Offset matrix A elements are changed as follows to apply to 
Channel 1 : 

Akn =0 for k=l 

30 - [0 0 ... V|/(-ML/2)r^ ... i|/( (L-k)M+l) W^^'^^] 

1 2 kM+1 LM+1 

which means the ACI error vector 5E is 
8E = L X 1 column vector 
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== A H W 



ACI metric for the two contiguous channels 
J(ISI) = 2 5E^ 8E Eigenvalue 
5 = Non-linear quadratic function of hf 

= 2 LS 



where the factor ""2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
10 3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 
allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI. 

15 

Cost function J for the LS algorithms is the weighted sum 
of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 
algorithms minimize J by selecting the optimal set of frequency 
coordinates {hf{k),Vk} for the selected set of parameters used 
20 to specify the characteristics of the dc waveform, frequency 
design coordinates, LS metrics, and weights. Cost function and 
optimization techniques are given by the equations 



25 Cost function J (27) 

J = ^^(/we/rzc^) J(metrics) 

X 

= weighted sum of the LS metrics J (metrics) 

where 

metrics = passband, stopband, deadband, ISI, ACI 
30 {w (metrics) , V/ne/r/c5'} = set of weights 

y^w(metrics) = 1 normalization 
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Optimization goal 

Goal: minimize J with respect to the selection of 
the {hf(k),Vk} 

5 

Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 
where the eigenvalue optimization algorithm uses the non- 
linear quadratic formulations of the LS metrics and the LS 
10 optimization algorithm uses the norm formulations for the 

LS metrics. 



FIG. '4 is a summary of the LS metrics and the construction 
15 of the cost function J. Design parameters 23 are the input and 
output design parameters. Input parameters are the number of 
polyphase channels M or equivalently the number of digital 
samples at spacing T over the symbol interval Ts=MT, the length 
of the FIR time response for the waveform in units of L which are 
20 the nunnber of digital samples per waveform repetition interval Ts 
so that the total number of digital samples for the symmetric FIR 
time response is equal to ML+1, number of DFT samples per FIR 
length n_fft for implementation of the LS algorithms, passband 
radian frequency <0p, stopband radiam frequency ©a, waveform 
25 repetition rate in radian frequency cOs, selection of the set of 
design coordinates {hf} to be used in the optimization, and the 
metric weights {w (metrics) } . Output parameters are the set of 
harmonic design coordinates {hf} that minimize J. Band metrics 

24 are the passband, stopband, and deadband metrics defined in 
30 equations (23), (23), (24) respectively. Interference metrics 

25 are the ISI and ACI metrics defined in equations (25) and 
(26) respectively. LS cost function J 26 is the weighted 
linear sum of the metrics defined for the band 24 and the 
interference 25 as defined in equation (27) . 
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FIG. 5 — i-9 — a — flow — diagram — — tfee — LS roouroivo — oolution 
algorithm. — There arc two loopo with the topology conotructcd oo 

that the outor loop — 27 io an iteration over the oot of metric 

5 weights — (w (mctrico) ) / — etftd — tfee — inner — noDted — loop 28- irS — tiie 

roouroivo or itcratived LS oolution to find the optimal — — 
the given deoign parametcro and weighto where each otep reforo to 

an iteration otep in tho oolution. A roouroivo LS oolution io 

required to be ablo to oolvo the highly non - linoar intcrforoncc 
10 motrico ao well ao the band metrico when one ohoooco to uoc non - 
linear tochniqueo to conotruot theoc metrioo. 

FIG. 5A,5B,. . .,SP,5Q are the Matlab 5.0 code for the LS 
recursive eigenvalue solution algorithm used to design the 

15 Wavelet Wavelet in FIG. 6. Included are a listing. of the frequency 
domain design harmonics, time response/ and data plotting 
software used in the optimization^ and an example of the scaling 
of this mother Wavelet into a Wavelet with new scaling 
(dilation), time translation, frequency translation parameters , 

20 and the functions used in the code. The Matlab code uses an 
iterative approach to find the optimal frequency domain 
eigenvector solution to the quadratic form for the LS error whose 
error matrix is the weighed sum of the error matrices for the 
stopband, passband, deadband metrics and the error matrices for 

25 the ISI and ACI metrics which are functions of the eigenvectors > 
For each iteration the ISI and ACI error matrices are constructed 
with the frequency domain eigenvectors from the previous 
iteration. The iteration finds the new frequency domain 
eigenvectors which are then used to update the ISI and ACI 

30 metrics for the next iteration. This iteration is repeated until 
there is convergence. The Matlab code allows the operator to 
update the metric weights as well as the scenario parameters to 
find the optimal choice for the performance goals. 
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FIG 5A Matlab code starts with the selection of the design 
parameters in 29 Step 1> Step 1.1 lists the scenario parameters 
which are M=16 digital samples per Wavelet \\/ sample interval, 
L-16 nominal length of Wavelet V[/ in units of N=ML+1 Wavelet \\f 
5 length which is represented by W in equations (10) , fs=l is the 
normalized channel spacing equal to the Wavelet \\f sample rate, 
fp=0.8864 is the passband frequency interval, n f-16 is the 
number of design frequency harmonics, n fft=1024 is the number of 
digital samples used for calculation of the frequency spectrum 

10 centered at 0 frequency. ebno=Eb/No==6. 0 dB is the design value of 
the energy per bit Eb to noise power density N p ratio, and 
X imbal aci=6.0 dB is the assumed channel-to-channel power 
imbalance used in the calculation of the ACI errors. Step 1.2 
derives the software parameters from the design parameters in 

15 Step 1.1. Step 1.3 lists the optimization parameters which are 
the number of iterations n iteration = 10, and the metric weights 
w (metric) from equation (27) and in FIG. 4 equal to w(pass) = 
alpha 1=1. e-2 for passband, w (stop) =alpha 2=0.80 for stopband, 
w(ISI)=alpha 3=2. e-3 for ISI, w(ACI)=alpha 4=0.5 for ACI, and 

20 w(dead)=alpha 5=0 for deadband. The number of iterations was 
selected to provide a convergent solution for the weighted LS 
error metrics w (metric) J (metric) in FIG. 5K 38 Step 10 figure (2) 
in the Matlab code and for their sum J in figure (1) in the Matlab 
code. Weight values were selected to optimize the Wavelet i|/ 

25 filter performance in FIG. 6 and in figure (3) in the Matlab code 
in FIG. 5L 39 Step 11, Wavelet \|; ripple, ISI, ACI signal-to- 
noise SNR power ratio losses in figure (4) in the Matlab code, and 
Wavelet \\f time response in figure (5) in the Matlab code. 

30 FIG. SB 30 Step 2 are the initialization calculations prior 

to the start of the iteration loop. Step 2.1 calculates the one- 
sided Wavelet y length parameter m=128. Step ,2.2 calculates the 
matrix transformation bw matrix which maps the \\f frequency domain 
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design eigenvectors into the one-sided v|/ time response. Step 2.3 
refers to the function pmn in FIG. SP^SQ 41 in Step 13.1 which 
computes the LS error matrix which is the sampled data integral 
Z[8c»8cnda) over the band in equation (23) for the the passband 
5 and stopband metrics J (pass) , J(stop) . The function pmn d 
computes the corresponding integral over the deadband in equation 
(24) and is not listed since it was not used for FIG. 6. Step 
2.4 constructs a structural matrix c matrix used in the ISI and 
ACI calculations. Step 2.5 constructs the sample rate templates 
10 for the passband upper edge, sample rate marker, and passband 
lower edge for frequency performance evaluation. 

FIG. 5C 31 Step 3 uses the function pmn to compute the LS 
error matrices passband and stopband in Steps 3.1, 3.2. Deadband 
15 LS error matrix deadband is set equal to a null matrix in Step 
3.3. Step 3.4 computes the matrix p total equal to the weighted 
sum of these matrices. Step 3.5 maps this LS error matrix 
p total in the time domain into the LS error matrix pw t in the 
frequency domain using the matrix transformation bw matrix. 

20 

FIG. 5C,5D 32 Step 4 initializes matrices used in the 
iteration loop and starts the iteration loop for the iterations 
i iteration=l,„,,10. Step 4.1 calculates the eigenvalue and 

25 eigenvector which minimize the cost function J in equations (27) 
and FIG. 4 whose LS error matrix is pw t. Step 4.2 uses the 
bw matrix to map the eigenvector in frequency into the Wavelet 
time response xy-hn and stores the eigenvector as the set of 
Wavelet V[/ design harmonics vt/k=hw eig(k) for k-0,1,..., 15 where h/ k 

30 has been defined in equations (9) and in equations (11) upon 
recognizing that there are n f=16=L design harmonics k o -k=0,l,. . 
.,15. There is no reason to use the negative harmonics since the 
spectrum is real and symmetric about the 0 frequency. Step 4.3 
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computes the passband peak-to-peak: ripple as xripple and the 
stopband attenuation as xstop. 

FIG. 5E 33 Step 5 calculates the weighted passband error 
metric w (pass) J (pass) =beta pass, weighted stopband error metric 
w(stop) J(stop)=beta stop, and weighted deadband error metric 
w(dead) J(dead) ==beta dead in equations (27) using the b vector 
format for hn and the LS error matrices passband, stopband, 
deadband. 



FIG. SE,5F 34 Step 6 calculates the ISI and ACI LS error 
matrices, their metrics, and their error contributations to the 
SNR loss. Step 6.1 calculates the ISI LS error matrix w matrix 
and the corresponding ISI metric J(ISI)=errM isi with Matlab code 

15 which implements equations (25) , and calculates the ISI residual 
error errV isi that causes an ISI SNR loss. Step 6.2 calculates 
the ACI LS error matrix w f matrix and the corresponding ACI 
metric J(ACI)^ errM aci with Matlab code which implements 
equations (26) , and calculates the ACI residual error errV aci 

20 that causes a ACI SNR loss. 

FIG. 5G 35 Step 7 calculates the weighted LS error metrics 
for ISI and ACI and the LS error matrix pw t for the next 
iteration. Step 7.1 calculates the weighted ISI LS error metric 
25 w(ISI) J(ISI)=beta isi and the weighed ACI LS error metric 
w(ACI) J(ACI) =beta aci in equations (27) and the sum of the 
weighted LS error metrics J=errM isi. Step 7.2 saves the weighted 
LS error metrics and their sum J for each iteration. Step 7.3 
updates the pw t matrix for the next iteration. 

30 

FIG. 5G,5H 36 Step 8 calculates the SNR losses, stores 
these losses for each iteration, and completes the iteration 
loop. Step 8.1 computes the SNR loss in dB due to passband 
ripple as xloss ripple, due to ISI as xloss isi, due to ACI as 
35 xloss aci, and the total SNR loss as xloss total. Step 8.2 saves 
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these SNR losses for each iteration and completes the iteration 
loop. 

FIG. 5H,5I,5J,5K 37 Step 9 documents the Wavelet \\f 
5 frequency design harmonics and time response generated by the 
Matlab code for the Wavelet in FIG. ^^6. Step 9.1 lists the 
Wavelet \\f frequency domain harmonic design coordinates 
\l/ k(k)=hw eig(k) for harmonics k=0, 1,..., 15. Only the positive 
harmonics are listed since the frequency spectrum is real and 
10 symmetric about k=0. Step 9.2 lists the Wavelet \|/ time response 
V|/(n)=hn(n) for n='0^1,2,. . . ,128 with the coordinate frame 
centered at the Wavelet peak. Only the positive digital samples 
are listed since the time response is real and symmetric about 
n=0. 

15 

FIG. 5K 38 Step 10 plots the total weighted LS error metric 
J in figured) and the individual weighted LS error metrics in 
figure (2) as functions of the iteration number. After about 6-7 
iterations both the total and the individual weighted LS error 
20 metrics are observed to stabilize indicating that the selected 

number of iterations n iteration=10 is sufficient - for 

convergence . 

FIG. 5L 39 Step 11 plots the Wavelet ij/ frequency response 
25 in figure (3) in the Matlab code which is the same as the 
frequency response in FIG. 6, plots the SNR losses due to 
passband ripple, ISI, ACI and the total loss in figure (4) in the 
Matlab code, and plots the time response in figure (5) in the 
Matlab code for use in the optimization of the LS metric weights 
30 in Step 1.3 as well as the scenario parameters in Step 1.1. 

FIG. 5M 40 Step 12 documents the Matlab code to rescale the 
mother Wavelet \\f in FIG. 6 which was derived in the previous 
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Matlab code in Steps 1-11/ for new scale (dilation) p, time 
translation q, and frequency translation k parameters compared to 
the values p=0,q=0,k=0 for the mother Wavelet y. An example set 
of new parameters documented in the code is p=2/ q^2, k=3 and 
5 with the assumption that the digital sampling rate remains fixed 
under the parameter change. A constant digital sample rate keeps 
the frequency band constant under the parameter change. 

FIG. 5M 40 Step 12.1 defines the new Wavelet v|/p , q, k in 
10 equations (7)^(18) resulting from these parameter changes, in 
equations (28) 

New Wavelet scaled from mother Wavelet (28) 

15 ilZr>.a.k = i|/(2'Pn-qM) exp(j27t k /ML) 

Case 1 Scaled (dilated) sampling with fixed M 

subsamples (decimates) n by the factor 2^ 

* Represent n as 



20 n = nO + 2Pnp 

n 0 = 0,1, . . 2P'l 

n p = 0, +/-lr +/-2, . . . 

' We observe 2"^n = n p 

' This enables the new Wavelet to be written 

25 M/o.a.k = vi/(n p-qM) exp(j2n k n p /ML) 

Case 2 Sampling n remains fixed and M is 
rescaled by the factor 2^ 



• New M is M new 



30 M new = 2^ M 



* This enables the new Wavelet to be written 
y[p,q,k = 2"^^^ i|/(n"qM new) exp(j27t kn/M newL) 
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wherein the new Wavelet equation from (7) ^ (18) is re-written 



using the frequency translation parameter k=r and the explicit 
frequency translation expression exp(j27t k 2"^n /ML) which is 
recognized as the discrete Fourier transform kernal. 

5 

In Case 1 in equations (28) corresponding to the standard 
Wavelet assumption, the subsampling or decimation of n by the 
factor 2^p stretches the time between contiguous samples by the 

10 factor 2^p while keeping constant the number of samples M per 
Wavelet v|/ sample interval, and therefore reduces the frequency 
band by the factor 2'^p, The new n can be represented as the sum 
of the two fields n 0 = 0,1,. > , (2^p -1) and n p = 0, +/-1, +/- 
2, . . . Scaling or equivalently subsampling results in the new 

15 n equal to n new = 2'^-p n = n p corresponding to a stretching 
of the sample interval by 2^p. The new Wavelet length N new 
remains constant N new=N using the n p sampling. 

In Case 2 in equations (28) which is the communications 
20 application being considered, the sampling is held constant 
corresponding to a constant frequency band assumption. The factor 
2''-p is divided out in the Wavelet expression which is equivalent 
to the M being scaled to M new=2^p M corresponding to increasing 
the number of samples per Wavelet interval by the factor 2^p. 
25 This means Case 2 corresponds to increasing the maximum number of 
channels by the factor 2^p over the same frequency band. The new 
Wavelet length is equal to N new=M newL+1 =1025 since 
M new=2^pM=64. 

30 FIG. 5M 40 Step 12.2 derives the new matrix transformation 

bw matrix new which . maps the \\f frequency domain design 
eigenvectors into the one-sided xjf time response. 
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FIG. 5M,5N 40 Step 12.3 uses the matrix transformation 
bw matrix new to map the design eigenvector hw eig into the new 
Wavelet in three steps. In the first step, the mother Wavelet is 
translated into the new baseband Wavelet Vi/D ,Q==o,k=o hn 0 which is 
5 the mother Wavelet scaled by p and with no time and frequency 
translations q=k=0. The second step translate this baseband 
Wavelet by qM new to generate vt/p ,q,k=o hn 1. The last step 
translates this Wavelet it/p ,q,k^o = hn 1 in frequency to give the new 
Wavelet M/p ,q,k = hn new. 

10 

FIG. 5N 40 Step 12.4 plots the Wavelet time responses for 
the mother Wavelet V|;=hn and the new Wavelet xt/p ,q,k = hn new in 
figure(6) in the Matlab code. FIG. 50 Step 12.5 plots the 
Wavelet frequency responses for the mother Wavelet i|;=hn and the 
15 new Wavelet vj/p ,q,k - hn new in figure (7) in the Matlab code versus 
the normalized frequency/ (\t/=hn sample rate) and in figure (8) in 
the Matlab code versus the normalized frequency/ (it/p ,q,k = hn new 
sample rate) . 

20 FIG. 5P 41 Step 13 documents the functions used in the 

Matlab code. Step 13,1 is the code for the function pmn which 
computes the matrix for the J(BAND) in equations (23). FIG. 5Q 41 
Step 13.2 lists the code for the function freq rsp which computes 
the Fourier transform of the input hn versus the normalized 

25 frequency/ (Wavelet v[/=hn sample rate). 

ConoidGr the innor nootod loop 2 8 . The rocuroivo LS - 
ee l-ution otarto with the initial otop i-^0 29 which bcgino with 
^c oclcction of the input dooign paromotoro and wcighto and tho 
30 ooloction of an initial oot of valuoo for the -j^ ^] 30. Next tho 
ha nd motrico arc calGulatcd 31. Since this ia the initial otop 
-3r ^ 32 the coot function J io rootrictcd to tho linoar band 
mctrico in order to fine tho linoar approximation to (h^ 33 to 



4-ft- itializo the non-linear oolution otarting with otcp i-1 3 4 
whorcin the highly non - lincar intcrforcnco mctrico arc included 
-ift J. Following the oignal flow/ for otep i-1 the band metrica 
33r and interforoncc mctrico 35 arc calculated/ weighted/ and 
5 oummod to form the coot function J 3 6 . An iterative oolut> eft 

algorithm 37 findo the boot approximation in otep i-1 to the 

•fi^d — which minimizeo J uoing the — approximation to — RirJ — from the 
prcviouo otcp i-0 — to linearize the search algorithm coefficiento 

for the eigenvalue and LS otep i-1 iteration. If there io no 

10 oonvcrgenoe — fee — the — correct — oolution — fof — [-fed — thio — rccuroivc 

oequence — — calculationo — is — repeated — §ef — otep — i=2 34-: Thio 

rccuroivc — oolution — technique — ts — repeated — — oubocquent — otepo 

until there io convergence — 38 — whereupon one exito thio inner 

loop. 

15 

Conoider the outer loop. After exiting the inner loop — 30 

fefee — solution — [it^ — — tooted — to — see — ar€ — it — meeto — the 

performance goalo 39-; If not, — a new oet of metric weights — i& 

oolccted and the — inner — loop — is — initialized i-^0 40 and the 

20 inner — loop — is — uoed — fee — find — the — next — oolution — §ej? — (h^) . Thio 

proccoo — continue o — until — tfee — performance — goalo — ete — met — 

adequated — approximated 39 whereupon the — algorithm io — exited 

with the — final oolution oet [h^ 44-^ 

25 Applications of this new invention to both communications 

and radar will be given using these example algorithms and other 
algorithms supported by this invention. These new waveforms are 
considered for the applications: l)to replace the square-root 
raised cosine waveform (sq-root rc) which is extensively used for 

30 the third generation (3G) CDMA communications, 2) to replace the 
Gaussian minimum shift keying (OMSK) waveform for constant 
amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 
aperture radar (RAR) applications. 

35 
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CDMA communications application for the current and the new 
3G CDMA considers a waveform designed with this new invention as 
a possible replacement for the sq-rt rc waveform with bandwidth 
expansion parameter a = 0.22 to a = 0.40. This notation means that 
5 for a = 0.22 the spectral efficiency is (symbol rate /bandwidth) = 
1/1+a =1/1,22 = 0.82 = 82%. A basic advantage of the waveform 
is the potential for a symbol rate increase within the same 
bandwidth with an increase in the spectral efficiency to « 100% 
depending on the application and operational constraints. The dc 

10 power spectral density or power spectriom (PSD) of the waveform is 
compared to the PSD for the sq-rt r-c in FIG. 6. Plotted are the 
measured PSD in dB 42 versus the frequency offset from dc 
expressed in units of the symbol rate 43 Plotted against the 
normalized frequency offset are the dc PSD for the new waveform 

15 44, the sq-rt r-c with a=0.22 45, and the sq-rt r-c with 
a==0.40 46. It is observed that the PSD for the new waveform 
rolls off faster that that for the sq-rt r-c which means that the 
new mr waveform will support an increased symbol rate for a given 
available frequency band while satisfying the inherent 

20 requirements for low ISI and MAI (multiple access interference) . 

Constant amplitude BEM application of the new waveform 
indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMSK. The GMSK finds 

25 applications for transmitters which operate their HPA(s) 
amplifiers in a saturation , mode in order to maximize their 
radiated power from the HPA{s), and which require a BEM PSD to 
avoid excessive spreading of the transmitted power. Simulation 
data for the dc PSD is plotted in FIG. 7 for the new waveform BEM 

30 and the GMSK. Plotted are the measured dc PSD in dB 47 versus 
the frequency offset from dc expressed in units of the bit rate 
48. Plotted against the normalized frequency offset are the dc 
PSD for the new waveform BEM 49 and the GMSK 50, for a length 
parameter L=10 where L is the length of the phase pulse in terms 
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of the phase pulse repetition rate. The significance of this 
example data is that the new waveform has the potential to be 
designed to offer a PSD which is less spread out than the current 
GMSK, and therefore an improved BEM waveform. 
5 \ 

Radar RAR and SAR application of the waveform indicates 
that it is a viable candidate to replace the current chirp 
waveforms for wideband signal transmission, when combined with 
pseudo-random phase codes. Results of the simulation for the new 

10 waveform and an unweighted frequency chirp waveform are given in 
FIG. 8. Plotted are the ambiguity function for the new waveform 
51 and the unweighted frequency chirp waveform 52, The dc 2- 
dimensional radar ambiguity function 53 is plotted as a function 
of the frequency offset in units of fTp and the time offset in 

15 units of t/Tc where Tp is the phase-coded radar pulse length or 
length of the phase code and Tc is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
Tg so that Tc = Tg. It is observed that the new waveform has the 
potential for significant improvements in the ambiguity function 

20 and by implication in the performance. 

Preferred embodiments in the previous description ie — are 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 

25 embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein fe«tand— is not to be accorded the wider 

30 scope consistent with the principles and novel features disclosed 
herein. 
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Currently amended Claims 
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APPLICATION NO. 09/829,118 

INVENTION: N|^ Multi-Resolution Waveforms 

INVENTORS: Urbain Alfred von der Embse 



CLAIMS 

WHAT IS CLAIMED IS: 

Claim 1. (currently amended) A method for designing mcano for 
^fefee — deoign — of now m ulti-resolution waveforms and filters , said 
method comprising: in tho Fourier domain with proportioa which 
means fpr generating complex Wavelet waveforms and filters, 
means for including a frequency translation property specified by 

a frequency translation parameter in addition to the 

scale and translation parameters, 

provide — cxtcnoiono — — the Wavelet — concept — to the — Fourier 

domain or oquivalontly the frequency domain 
means for using a subset of the Fourier harmonics as the design 

coordinates (harmonics) specifying the waveform design, 
means for p rovide p roviding a single waveform designs- for all 

of the w avefor ms at multiple scales, translations, and 

frequencies, and o at multiple ocaloo 
means for using the design harmonics and frequency translation 

property to generate Wavelet waveforms at multiple scales 

and multiple frequencies. 

provide — dcoigno — which — rcquirco — a — relatively — few — dcoign 
coordinatco — compared to — fefee — FIR time — rooponoo — oamploo — 
the waveform 

provide design — methodologies which een incorporate 

application mctrico to improve tho waveform performance 
provide — dcoigno — that — allow — the — use — of — direct — dcoign 
mcthodologico — that — circumvent the need to — solve a Wavelet 
Iterated filter bank construction to obtain the wavoformo 
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thoroby providing improvod flexibility fee moot the 

application goalo 

Claim 2, (currently amended) g -= — A mcanp — for the dooign of multi 
rcoolution — wavoformo — aftd — f iltcro — ift — fcfee — frequency — domain — ^ 
equivalontly the Fourier of the t - f opaco for waveform and filter 
gpplicationo, with propertieo which The method of claim 1 further 
comprising: 

means for combining filter design requirements with Wavelet 

requirements to generate finite impulse response filters 
(FIR) and waveforms ^ 

provide a meano for oingle deoign to be uoed for all of the 
wavoformo at — the multiple reoolutiono or frequency bando 

provide a meano — for dooigno — that require a relatively few 

dooign — coordinateo — compared — fee — fefee — number — ef — digital 

oamplco covered by the Wavelet 
means p rovide a meano for the deoigno to for include including 

frequency and 

time application — error m etrics fee — improve — fei=5e — waveform 
performance and for including application 

constraints in the weighted least-squares sum J of the 
error metrics, 

means for finding the FIR design harmonics that minimize the 

constrained sum J of the error metrics , 
means for implementing direct least-squares error techniques to 

find the design harmonics that minimize J, 
means for implementing iterative eigenvalue error techniques to 

find the design harmonics that minimize J, and 
means for using the plurality of design methodologies in the 

time- frequency space to find the design harmonics. 

provide — a meano — to uoc direct dooign methodologioo — in the 

frequency - time — domain — c for — aH — — fefee — wavoformo — aftd 

filters for the multiple frequency bando 



2 



Claim 3. (currently amended) The methods of claim 1 and claim 2 
further comprising; 

means for designing polyphase orthogonal filter banks with 

almost perfect reconstruction properties, 
means for designing polyphase quadrature mirror filter 

(QMF) orthogonal filter banks, 
means for designing waveforms and filters with no excess 
bandwidth, and 

means for designing almost ideal waveforms and filters for linear 
and non-linear applications. 

A method for tho dcaign of multi - rooolution wavoformo which 

allows — Fourior — domain — tochniquoo — to bo — used. Proportioo 

eeft — inoludo — oomo or moro — of tho — liotcd proportioo — iB — 
and (2) 



Claim 4. (currently amended) The methods of claim 1, claim 2, 
and claim 3 further comprising: 

means for including analytical, iterated filter bank, and 

scaling function design techniques for tiling (covering) , 

means for tiling the time-frequency space for 
communications and radar, 

means for tiling the spatial-frequency space for media 
processing, 

means for tiling the time-frequency-beam space for cellular and 

satellite communications and radar, and 
means for tiling the time-wavelength space for laser and optical 

communications, targeting, and ranging. 

A method for the dcoign of multi - rooolution wavoformo which 
eaR — incorporato — Fourior — domain — toohniquco — into — dooign 
mcthodologioo — which — eaR — include — analytical — aftd — itoratod 
filter bank oonotruction design tcchniquco. 
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Claim 5. (currently amended) A method for the design of multi- 
resolution waveforms and filters^ said method comprising: 
means for generating complex waveforms with frequency 

translation^ scaling, and time translation properties^ 
means for providing a single waveform design for all waveforms 

at multiple scales, translations, and frequencies, 
means for using a subset of Fourier harmonics as the design 

coordinates specifying the waveform design, 
means for using the design harmonics and the frequency 

translation property to generate waveforms at multiple 
scales and multiple frequencies, and 
means for using the plurality of optimum design techniques to 

derive the design harmonics. A method for tho analyoio and 
dcoign of multi - rooolution wavcformo uoing Fourier domain 
tQchniquco which take advantage of the new invention 
diocloourco on tho charactoriaation and dcoign of multi - 
rcoolution wavoformo in tho — Fourier domain. 



Claim 6. (currently amended) The methods of claim 5 further 
comprising: 

means for designing the multi-resolution waveforms for bandwidth 

efficient modulation (BEM) , 
means for designing the multi-resolution waveforms for high power 

amplifiers and non-linear signal processing, 
means for designing the waveforms for non-linear applications, 
means for designing the waveforms for multi-resolution and 

bandwidth-on-demand communications , . 
means for designing the multi-resolution waveforms for synthetic 

and real aperture radar applications, and 
means for designing the multi-resolution waveforms for laser 

radar and laser communications. A now formulation in 

-(-&-) — for multi - rooolution waveform ao a function of tho dc 
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multi - rcoolution — wavoform — which — addo — the — oonccpt — — a- 
frequoncy indox — that allowo — the multi-rooolution waveform 
to be placed arbitrarily throughout the t^f opaoo thoroby 
ir) — avoiding the rcotrictiono of the Wavelet iterated filter 
conotruction — fe*^ — tiling — a — # — opaco/ — aad — 2-) — allowing — fefee 
Bew — multi - rcooloution — wavof ormo — fee — be — uood — So^ — multi - 

reoolution — coinmunicationo aftd for bandwidth on ■ demand 

communicationo application, in place traditional 

Wave Ic to . 
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FIG. 3 PSD Requirements for Coimnunications 
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FIG. 4 LS tfetrics and Cost Function 
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Design Parameters 



Passband Eq. 23 
Stopband Eq. 23 
Deadband Eq. 24 



ISI Eq. 25 
ACI Eq. 26 



LS Cost equation 
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FIG. 5 LS Recursive Solution Algorithm 
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Nozmalized freqaencs', syadsol rate 43 



5 



10 



15 



54 



FIG. 7 PSD for New Waveform BEM and GMSK 
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FIG. 8 Radar AoibiguitY Functions of New Waveform 
and Unweighted Chirp Waveform 



New Wavelet Ambiguity FimctiOT 51 Unweigjited Chirp Ambiguity Functi(Mi 52 
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